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The translational diffusion coefficients of the 45S and 50S states of the large ribosomal subunit o£E. coli were 
determined from the spectral distributions of quasielastically scattered light with a 5 mW He-Ne laser as the 
source. The spectral analysis was performed l>\ directly Fourier transforming the photocurrent and fitting to 
double-Lorentzian profile via a non-linear regression routine. A small amount (about 1%) of strongly scattering 
contaminants required the double-Lorentzian profile in order to extract the diffusion coefficients of the principal 
components. The results are: D 20 °, w(45S) = (1.79 ± 0.12) X 10" 7 cm 2 /sec, and D 2() °, w(5()S) = (1.91 ± 0.06) 
X ICT'enr/sec, the latter being in accord with those reported in the literature. The transition from the SOS state 
to the 45S state is not attended by a change in its molecular weight il the partial specific volumes are assumed to 
be the same. 
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1 . Introduction 

The most immediate task in characterizing complex sub- 
cellular particles such as ribosomes, is to assign molecular 
weights. This is commonly performed by determining; the 
sedimentation and diffusion coefficients in conjunction with 
the partial specifie volume. With the advent of the laser light 
scattering technique, [l] 1 the diffusion coefficients of a 
variety of subcellular particles, viruses, and macromolecules 
[2] have been determined at a substantial savings in effort 
and time, though not necessarily with improved precision 
over conventional techniques. As diffusion coefficients be- 
come more readily available, molecular weight assignments 
will be less time consuming. Interest in the diffusion coeffi- 
cient also resides in deducing the geometric eccentricity [3] 
of a uniform density particle by comparing the equivalent 
Stokes radius to the radius of gyration obtained by the 
scattering intensity technique with X-ray, slow (or cold) 
neutron, or visible light. Conversely the component density 
segregation within a particle can be inferred when the Stokes 
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radius deduced from the diffusion coefficient is larger than 
the radius of gyration. Such an inference was drawn by 
Serdyuk el al. [4] with the SOS subunit from the study of 
small angle X-ray scattering and the diffusion coefficient. 

In this paper, we report a study of the diffusion coefficients 
of the 50S and 45S states of the large ribosomal subunit 
using quasielastic light scattering (QEL). The purpose is two 
fold. First, we show that the two states, which differ 1>\ 5 
Svedberg units, can be distinguished by the diffusion coeffi- 
cient as obtained by the technique. Second, a method for 
extracting individual spectral halfwidths, hence diffusion 
coefficients, from a bimodally distributed scattering system 
is put forth. The scheme is analogous to the one adopted b\ 
Bargeron [5] except that ours is in the frequency domain 
while his is in the time domain. 

In the determination of the translational diffusion coeffi- 
cient by QEL, a trace amount of contaminants having a 
larger dimension than that of the principal component poses 
a major experimental problem. This is simply due to the 
larger scattering intensity of the impurity component which 
obscures the scattering contribution of the principal one. 
This sort of polydispersity effect has been the cause of a good 
deal of examination in test systems recently [5-8]. It is 
generally agreed that one indeed measures the z-average 
quantity [8-10] of the diffusion coefficient when the intrapar- 
ticle form factor is unity as first shown by Koppel [11]. 
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Hence the data analysis of the spectral profile of scattered 
light from a bimodally distributed system compared to a 
smoothly distributed system requires some care, particularly 
when one is interested in a small component that is contam- 
inated by even trace amounts of a large component. Depend- 
ing on the ratio of linear dimensions of the two, one can 
easily be misled by the z-average quantity that is predomi- 
nantly weighted by the larger component [12]. We show here 
that the spectral analysis has been successfully carried out 
in determining the diffusion coefficients of the 45S and 50S 
states of the large ribosomal subunit contaminated by some 
larger particles, from sources as yet undetermined. The 
contaminant components appeared to be 6 to 7 times greater 
in diffusion coefficients than the principal ones, but present 
in no more than a few parts per hundred by weight. 



2. Experimental 

2.1 Sample preparation 

The ribosomal subunits were extracted from the MRE-500 
strain of E. coli in Professor Spirin's laboratory at the 
Institute of Protein Research, The Academy of Sciences, 
U.S.S.R. by a modified version of the method reported 
previously by Gavrilova et al. [12] with the use of B-XV 
zonal rotor of a MSE* superspeed 65 preparative ultracentri- 
fuge. All solutions were prepared immediately before the 
experiments from the ribosome preparations stored in satu- 
rated ammonium sulfate. The 50S ribosomal subunit solution 
was prepared as follows. The ribosome suspension was 
centrifuged for 25 minutes at 1.7 X 10 4 g to pellet. After the 
supernatant liquid was decanted, the pellet was resuspended 
in a few ml of the buffer (0.07 M KC1, 0.01 M tris-HCl, 1 
mM Mg +2 , pH 7.1) and dialyzed for eight hours against the 
same buffer with threefold exchange. The 45S state of the 
subunit was prepared from a 50S subunit solution by dialysis 
against the following buffer: (0.5 M NH 4 C1, 0.01 M tris-HCl, 
1 mM Mg +2 , P H7.1). 



2.3 Light Scattering 

The instrument for this study has been described else- 
where [13]. The laser source of the instrument is a 5 mW 
He-Ne laser (Spectra-Physics 135)*. The time varying output 
of a photomultiplier (RCA C7164R)* produced by the 
homodyne beating of the scattered optical field at the 
photocathode [14] is converted to voltage, amplified by a 
wide-band device, band limited to frequencies below the 
Nyquist limit of one-half the sampling frequency by a filter 
(Dytronics 722)*, sampled and digitized at an appropriate 
rate controlled by a real time clock, then directly Fourier 
transformed, and accumulated on an on-line minicomputer 
(PDP8/E)*. Over a given bandwidth (2.5 KHz and 10 KHz 
for this study), the power spectrum is constituted of 256 
power coefficients that are averaged over 2000-4000 scans. 
After correcting for the instrumental profile, the spectrum is 
then analyzed by fitting it to model profiles via a non-linear 
regression routine [15]. 

The scattering angles were incremented by 15° from 45 — 
90° and finally by 20 to 110° for both 45S and 50S solutions, 
and an additional angle at 30° for the latter was examined. 
Temperature of the scattering cell (10 mm NMR tube) 
immersed in a refractive index matching liquid (Silicone oil) 
was not controlled away from the ambient but continually 
monitored with a calibrated quartz thermometer (Hewlett- 
Packard 2801 A)*. For a given run of scattering, the temper- 
ature variation was noted to be within 0.2 °C. 



2.4 Data Analysis 

For a bimodally distributed scattering system, the scat- 
tered power spectrum in the frequency domain is given by 
the following equations, [14] provided (1) the optical field is 
governed by Gaussian statistics and (2) the detection scheme 
is through homodyne beating. 
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2.2 Sample Clarification 

The subunit solutions at a nominal concentration of 2 
mg*ml -1 were centrifuged on a preparative centrifuge for 1— 
2 h at 6.7 X 10 3 g before subjecting them to light scattering. 
In one case, we further clarified the solution by filtering 
through 0. 1 jLtm Nuclearpore* filter without imposing a 
positive pressure gradient upstream. Although the Schlieren 
patterns of the solutions before and after the scattering 
invariably produced a homogeneous single peak, the spectral 
analysis required the double-Lorentzian scheme (see below). 



where the current autocorrelation function is expressed as 



C t (T) = e < i > 8(t) + < i 2 > g {2) (r) 



e\r) 

? u, (r) 



1 + |g">(T)| 2 

/ ie - r ' lTl + / 2 e-^ lTt 
r, = D t q\ q = (4ir/X')sin(0/2) 



(2) 
(3) 
(4) 



Here, the electronic charge is denoted by e, the mean 
photocurrent < i >, Delta function ${t), the mean square 



photocurrent < i 2 >, and the normalized intensity autocor- 
relation function g (2) (r) is given by the Siegert relation, [16] 
eq (3). The normalized field autocorrelation function # (1) (t) 
is expressed by a linear combination of two exponential 
decays with the time constants T x and T 2 , which are 
proportional to the respective diffusion coefficients, D x and 
D 2 , and the corresponding intensity factors I\ and 7 2 . The 
wavelength of the light in the scattering medium is denoted 
by A. . Combining eqs. (1) through (4), we obtain the power 
spectrum in the broadened frequency scale from the incident 
angular frequency Cl) , i.e., v = (a> — (O )/27T; 

P t (v > 0) = e < i >/tt + < i 2 > 8'(v) 
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where 8'(^) is defined through 

b'(v)dv = 1 



The first term in the RHS of eq. (5) represents the shot-noise 
level, the second term the dc contribution, and the third the 
double-Lorentzian spectral profile. Thus fitting to the 
double-Lorentzian profile of an observed power spectrum 
involves five parameters; they are (1) the shot-noise level, (2) 
the constant total intensity factor of the spectral portion 
I 2 < i 2 >/tt 2 , (3) the halfwidth of the first component 
A^i/2 , (4) the halfwidth of the second one A^ 2 , and (5) the 
intensity ratio of the two components lijl\. We adopt the 
convention for designating the halfwidth of the slower moving 
component as Api/£ and of the faster one as Av^l so that 

A*81 > MK. 

The calibration of the double-Lorentzian analysis scheme 
was accomplished as follows. A binary mixture of polystyrene 
latex spheres was prepared by combining spheres having 
nominal diameters of 0.091 /xm and 0.810 fxm, in weight 
ratio of 98/2, i.e., 2 percent of the large spheres contaminat- 
ing the principal component. The experimental power spectra 
acquired at different angles were analyzed according to eq 
(5) and A^i/2 and APi/£ were determined. Four sets of the 
halfwidth values obtained at four different scattering angles 
provided two linear plots from which the diffusion coefficients 



of the two spheres were determined according to eq (6), and 
found to be in complete agreement with those predicted by 
the Stokes law within 3 percent. A parallel experiment with 
16 percent of the large sphere gave the analysis result which 
was indistinguishable from a solution consisting of pure 
0.819 jLtm latex spheres; the single Lorentzian analysis was 
good enough. Thus we can extract the halfwidth of the 
smaller component of a binary mixture if the larger particles 
are present in about 2 percent or less when their diameter 
ratio is about 1/9. 



3. Results and Discussion 

The spectra at 45° of the 45S and 50S preparations are 
displayed in figures 1 and 2 respectively. In each figure, the 
upper part shows the experimental data overlaid on the fitted 
single-Lorentzian (chained curve) and the double-Lorentzian 
(solid curve) profiles, and the lower part indicates the 
corresponding normalized residues for the two types of 
profile. The residue plots clearly demonstrate the successful 
representation of the observed data by the double-Lorentzian 
scheme in each case. Similar analysis was effected for all 
observed spectra and the same conclusion was drawn for 
every spectrum. In tables I and II we collect the halfwidth 
data together with the intensity ratios (/2//1) extracted from 
the double-Lorentzian analysis. Upon reducing all halfwidths 
to those in water at 20 °C, we determined the diffusion 
coefficients of the 45S and 50S states of the large ribosomal 
subunit and those of the contaminant particles via eq (6) at 
the hydrodynamic standard condition. In figures 3 and 4 we 
plot the reduced halfwidth of the principal components at 
different angles against r/ 2 for 45S and SOS preparations 
respectively. The error bar at each scattering angle repre- 
sents the ranges of A^i /2 required for the 95 percent 
confidence limit of the fitted spectral profile. In addition, we 
also present in Fig. 4 the corresponding plot for the contam- 
inating component of the SOS snbpartieles; the error bars are 
not shown for this case because they arc comparable to the 
size of halfwidth data points (asterisk). From these plots we 
calculate the diffusion coefficients of all four components: 

#20°,tt> (45S) = (1 .70 ± 0. 12) X MP 7 enr/sec. 
D 20 o 3W (45S cont.) = (2.22 ± 0.30) X K)" 8 cm 2 /sec. 

#20°,m> (50S) = (1.91 ± 0.06) X 10" 7 cm 2 /sec. 

D 2 o°,w (50S cont.) = (2.99 ± 0.14) X -10~ 8 cm 2 /sec. 

The uncertainty in each case is estimated by the 95 percent 
confidence interval for the mean value from the estimated 
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FIGURE 1. Upper: The observed power coefficients at 45° scattering angle of the 45S 
preparation vs the broadened frequency v(Hz); solid and dashed curves represent the fitted 
double -Lore ntzian and single -Lore ntzian profiles respectively. 

The ordinate is in arbitrary units. 

Lower: The normalized residuals for the two analysis schemes vs v(Hz); X for the double- 
Lore ntzian and + for the single -Lore ntzian scheme. 



standard deviation and Student's t- value at a = 0.05. The 
concentration dependence of the diffusion coefficient was not 
examined since its absence up to 0.6 mg-ml -1 for the subunit 
has been verified by Koppe [17] and it seems reasonable to 
assume that similar behavior would apply to the 45S state. 

Setting aside for the moment the nature and sources of the 
contaminant particles, we proceed with the discussion of the 
ribosomal subunits. The mean value of D2o°,u? f° r ine 50S 
subunit is in good accord with those obtained by Tissieres et 
at. [ 18] and Serdyuk et al. [3] with use of the conventional 
techniques, as well as that by Koppel [17] who used the same 
technique as ours. On the other hand the uncertainty in our 
determination is greater than the one in Koppel's work by a 
factor of two; we estimate 3.2 percent whereas Koppel quotes 
1.6 percent. The difference is attributable to the purity of his 
simple which, unlike ours, appeared to be free of contami- 



nants; he observed single exponential decay of the photocur- 
rent autocorrelation function which is equivalent to single 
Lorentzian profile in frequency domain. Thus our study 
seemingly attests the success of Koppel's ingenious scheme 
to effect the scattering with a pure sample at varying 
concentrations within a single cell. It is interesting to note 
that the two different signal processing routes, direct Fourier 
transformation in our case and photon count autocorrelation 
in Koppel's case, gave rise to an identical result although we 
could perform detailed statistical analysis with 200 data 
points for each spectrum while similar analysis might not 
have been possible with 20 data points (in a 24-channel 
correlator) by Koppel. In terms of the signal to noise ratio, 
the two procedures have been proven to be equivalent [19]. 

To the best of our knowledge, the diffusion coefficient of 
the 45S state of the large subunit has not been reported 
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Figure 2. Upper: The observed power coefficients at 45° scattering angle of the 50S 
preparation vs the broadened frequency v(Hz); solid and dashed curves represent the fitted 
double-Lorentzian and single-Lorentzian profiles respectively. 

The ordinate is in arbitrary units. 

Lower: The normalized residuals for the two analysis schemes vs v( Hz); X for the double- 
Lorentzian and + for the single-Lorentzian scheme. 



heretofore in the literature. The uncertainty in the final value 
however seems unusually large, amounting to nearly 8 
percent. Sources of the large error bound are unclear; the 
intensity ratios in the 45S are quite comparable to those in 
the 50S subunit which would rule out the possibility of 
greater contamination in the 45S case. A question we now 
face is: Is the difference in the mean values of the two 
coefficients [D(50S)-D(45S) = 0.12 X 1(T 7 cm 2 /sec) statis- 
tically significant? To this end we evaluate a 95 percent 
confidence interval for the true underlying difference [20]. 
We find (0.013, 0.227), which does not include zero. 
Therefore, we may conclude that difference between the 50S 
and 45S means is positive, i.e., that the difference is 
statistically significant at the 5 percent level of significance. 
Accepting the difference to be significant, we might conclude 
that the molecular weights of two states are the same because 



of the same ratio of sedimentation to diffusion coefficient 
provided that the partial specific volumes do not differ. Hill, 
Anderegg and van Holde [2.1] however found that the two 
states had different partial specific volumes but had the same 
molecular weight. Our results thus appear to contradict one 
of these two findings; either molecular weighl is different or 
partial specific volume is identical. A closer analysis how- 
ever yields consistency. Based on the error bounds of our 
diffusion coefficients and those of their partial specific 
volumes, we calculate according to the Svedberg equation 
the molecular weights of two states as (1.73 ± 0.16) X 10 (> g/ 
mol (5 2 o,m> = 45.2 X 10~ 13 sec) and (1.55 ± 0.07) X 10" g/ 
mol {s 2 o,w = 50.2 X 10~ 13 sec) respectively for 45S and 50S 
subunits. Hence, the two molecular weights are the same 
within experimental error. We should note here that the 
molecular weights, by sedimentation equilibrium, reported 



Table 1. The results of double-Lorentzian analyses for the 45S state of the large ribosomal subunit solutu 



Angle 
(degree) 


Temp 
(°C) 


(Hz) 


(Hz) 


to<$l/Av$l 


hlh 


(a) 

(/yPi) 8 


(b) 

(^V^Ocalc 


45° 


27.3 
(20.0,w) 


85.0 

(72) 


736 
(621) 


8.63 


0.397 


1.21 


1.21 


60° 


27.3 


159 
(133) 


1150 
(971) 


7.30 


.568 


1.73 


1.39 


75° 


27.1 


231 

(195) 


1754 
(1487) 


7.63 


.745 


2.27 


1.64 


90° 


26.8 


329 
(281) 


2452 

(2093) 


7.45 


.882 


2.69 


1.96 


110° 


26.3 


340 

(293) 


2843 

(2453) 


8.37 


.980 


2.99 


2.51 



After 0.1 jxm Nuclearpore'f filtration: 

60° 25.8 165 

(144) 



1271 
(1109) 



7.70 



1.400 



(a) Scaled (/ 2 //,) to match (P 2 /^i)caic at 45' 

(b) Calculated from P(q, R) = -j-$ 



(sin qR — qR cos qR) 



for a spherical particle with a radius R and the Stokes law, R — 



kT 
67717 qD 



where kT has the 



usual meaning, 170 is the viscosity of solvent, and/) the observed diffusion coefficient. 



TABLE 2. The results of double-Lorentzian analyses for the 50S state of the large ribosomal subunit solutu 



Angle 
(degree) 


Temp 

(°C) 


A^V/2 

(Hz) 


(Hz) 


A?»/MS 


hlh 


(a) 


(b) 
(P 2 /Pl)calc 


30° 


23.6 

(20.0,w) 


33.2 
(30.5) 


313 

(287) 


9.41 


0.313 


1.05 


1.05 


45° 


23.8 


86.5 

(79.0) 


679 
(620) 


7.85 


.528 


1.77 


1.11 


60° 


24.3 


184 
(166) 


1180 
(1066) 


6.41 


.669 


2.24 


1.19 


75° 


24.7 


278 
(249) 


1695 
(1517) 


6.10 


.836 


2.80 


1.30 


90° 


25.3 


403 
(355) 


2536 

(2239) 


6.30 


.873 


2.93 


1.43 


110° 


25.7 


504 

(440) 


3219 
(2816) 


6.39 


.902 


3.03 


1.63 



(a) Scaled {I 2 /h) to match (P 2 /Pi)calc at 30° 

(b) As in table I. 



by Hill et al. are (1.58 ± 0.08) X 10 6 g/mol (45S) and (1.55 
± 0.05) X 10 (i g/mol (50S), whereby the absolute values of 
molecular weight by two experiments are also in agreement. 
From these considerations, it is quite plausible that the two 
states have the same mass but differ in their conformations 
such that a slight change in geometric eccentricity (or axial 
ratio) is reflected in their hydrodynamic parameters. 

We now turn to the contaminants. As indicated above we 
really do not know where these come from. They might be 
the unfolded forms or large aggregates of the subunits though 
so small in amount that the Schlieren optics was not sensitive 
enough to detect them. The fact that the impurity component 
was representable by a Lorentzian profile (with A^i 1 /^) cannot 
constitute the necessary condition for its monodispersity. 



Any collection of components if present in small amount 
should have been fitted by a seemingly well defined one in 
the double-Lorentzian scheme, particularly when their diffu- 
sion coefficients are far different from that of the majority of 
scattering elements. Those solutions filtered through 0.1 jum 
Nuclearpore* filter gave rise to a larger intensity ratio /2//1, 
indicating that the contaminants could be removed, at least 
partially, upon filtration; they were not in equilibrium with 
the principal component. An example of this is provided in 
table I. While we cannot be sure of its character or sources 
without further study, it is possible to estimate its approxi- 
mate mass ratio to the principal component. In addition, we 
can say something about its axial ratio compared to the 
subunit particles. 
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FIGURE 3. The reduced spectral halfwidth (20°C, water) for the 

principal component of the 45S preparation is plotted against </~/tt. 
The error bars cover the 959£ confidence interval about the mean halfwidth determined bj ilu- 
double-Lorentzian method at each scattering angle. 




Figure 4. The reduced spectral halfwidth for the principal and 
contaminant components of the 50S preparation are plotted against q' 2 /7T. 

The error bars in the Av[% plot cover the 95% confidence interval about the mean halfwidth 
determined by the double-Lorentzian method at each scattering angle. 



weight, Ci the concentration in mass per unit volume, P{ the 
particle scattering function and Ri the linear dimension of 
the same component. In order to estimate the mass ratio 
Ci/c 2 , we need to make the following assumptions. They are: 
(1) K 2 = K u (2) Pi(q, R t )~ 1 at 30° ~ 45° scattering angle 
region, which might he low enough to allow such an 
assumption, and (3) M 2 is proportional to /)j~ 3 , which is 
equivalent to adopting the solid sphere model. From the last 
assumption, M 2 /M x arc estimated as 1.9 X ]()~ :i and 3.8 X 
10~ for the 45S and 50S states, respectively. The mass 
ratios c x /c 2 are then computed from e<] (7) as 5.8 X 10 -3 
(-0.6%) and 1.3 X K)" 2 (-1.5%) for the 45S and SOS 
preparations, respectively, indicating only a small fraction of 
the large contaminants in each preparation. A similar situa- 
tion was encountered by Chen el al. [23] in a recent study of 
dipalmatoyl phosphatidylcholine vesicles. 

The relative shapes of the principal and contaminant 
components may be deduced as follows. Since the only 
scattering angle dependent term in eq (7) is Pi(q, Ri), the 
observed / 2 //i ratio should be directly proportional U)P 2 /P\. 
Meantime, one can calculate the ratios P 2 /P\ assuming t lie 
sphere model for both components. We then compare them 
to the intensity ratio I 2 /l\; we scale l 2 /l\ so as to match the 
calculated P 2 /P\ at the lowest angle and examine the angular 
dependences of the scaled ratio, {P 2 /P 2 ) s vis-a-vis the 
calculated one, ( P 2 /Pi) Vi{ \v These are collected in the 
tables. It is apparent that (P2/F* i)s increases with the 
scattering angle faster than (^/^Ocalc? particularly with the 
SOS preparation. Thus we conclude that the principal ribo- 
somal subunit component is more asymmetric than the 
contaminant component because the particle scattering factor 
P(q, Ri) for a sphere decays most rapidly with q among the 
uniform density particles [24) and any asymmetry in the 
particle shape necessarily reduces the angular dependence 
of P(q, R 2 ). Our conclusion however must be moderated if 
the two components were to differ in the uniformity of particle 
density. 
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K 2 M 2 c 2 P 2 (q, R 2 ) 



(7) 



where K 2 is proportional to the refractive index increment of 
the ith component (over the buffer), Mi is the molecular 
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